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C ' Abstract Evidence for an accelerated expansion of the universe as it has been re- 

^ vealed ten years ago by the Hubble diagram of distant type la supernovae represents 

one of the major modern revolutions for fundamental physics and cosmology. It is yet 
unclear whether the explanation of the fact that gravity becomes repulsive on large 
scales should be found within general relativity or within a new theory of gravitation. 
However, existing evidences for this acceleration all come from astrophysical observa- 
tions. Before accepting a drastic revision of fundamental physics, it is interesting to 
r N ' critically examine the present situation of the astrophysical observations and the pos- 

• , sible limitation in their interpretation. In this review, the main various observational 

1-5 ' probes are presented as well as the framework to interpret them with special attention 

to the complex astrophysics and theoretical hypotheses that may limit actual evidences 



Q ■ for the acceleration of the expansion. Even when scrutinized with sceptical eyes, the 

evidence for an accelerating universe is robust. Investigation of its very origin appears 
as the most fascinating challenge of modern physics. 
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C^^ ' Modern cosmology has achieved remarkable progresses during the last fifty years. The 

C^ I general picture originally designed as the "Primeval atom" by Lemaitre and which 

\/-^ • has become the "Big Bang" model according to the word of one of its most famous 

c7^ ' opponent, F. Hoyle, is now recognized as the successful scientific representation of the 

^^ . world at the large scales (in space and in time) we can measurqj. Construction of 
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^ It is fair to say that few scientists are still opposed to the "Big Bang" picture. Most of the 
serious opponents try to demonstrate that some observational facts, most often only one, which 
are coming in support of the Big Bang may be interpreted in a different way and therefore 



this picture has necessitated the successive abandonment of philosophical and scientifi 
ideas, some of which are not only those of physicists but are also shared by the more 
general public. Maybe the first to be given up was the idea that it is hopeless to try 
to figure out a global picture of the universe. Although heroic pioneers can be traced 
back since a long time ago, undoubtedly Einstein is the first one to directly address the 
question of handling the universe as a single physical object with an appropriate tool 
in hands, general relativity (GR). The fact that one of the very first applications of his 
theory was Cosmology is an evidence that the cosmological question was central in his 
thoughts. This is reinforced by the fact that he proposed to modify the initial formula- 
tion of his theory given the problems he encountered. It is commonly said that Einstein 
introduced the cosmological constant in order to obtain a static universe because he 
was reluctant to the idea of an expanding picture. This formulation is inappropriate: 
it merely suggests that Einstein was willing to avoid an expanding universe, while he 
actually wanted to find at least one solution to the cosmological problem with an initial 
formulation which assumes the universe to be stationary, there is no indication in his 
seminal paper that he wished to reject a non-stationnary universe: he actually started 
his paper by discussing the problem of having mass at large distances in the Newtonian 
approach, noticing that this leads to divergence of the potential. He insisted that this 
would lead to unacceptably large velocities for stars. He also quoted that this can be 
cured by assuming a correcting term to the Newtonian potential equation: 

V^V5 - A((5 = 47rGp (1) 

and then proposed to modify his initial theory with the addition of the cosmologi- 
cal term A. This allowed him to construct the first relativistic cosmological model, 
the Einstein solution, which is spatially closed (because being spherical) and static. 
In 1919, de Sitter discovered a new solution to Einstein equation which was written 
in a stationary forrrjj and contains no matter (but a non-zero cosmological constant). 
It is only a few years after, in 1925, th at L emaitre identified the de Sitter solution 
to an homogeneous expanding universe l|l40l ). Friedmann (|9ll : |93) found the general 
homogeneous solutions, providing the equation for the scale factor R{t) and recognised 
their expanding nature. It is somewhat surprising that his work has remained totally 
unnoticed, despite a controversy with Einstein. During this period it is clear that the 
nature of the redshift of what Hubble had identified as extragalactic nebulae became a 
question addressed by many astronomers. Slipher's discovery was probably much more 
intriguing now that the nature of nebulae had been identified. Eddington is often men- 
tioned as the first astrono mer to have noticed a possible connection between redshift 
and the de Sitter solut ion (Il40 ''l . Carl Wirtz was clearly looking for a relation between 
redshift and di stan ce (|238r) and had in mind the possible cosmological information it 
might provide l|239r) . Lemaitre reestablished the equations Friedmann derived, showing 
that expanding solutions were leading to a redshift proportional to the distance. He 



the whole construction has to be questioned. It is useful to remember that the success of a 
scientific model is - in some sense — measured by the number of predictions it leads to and 
how many are successful. Newton theory of gravity is wrong, but nevertheless it remains a 
high quality scientific theory because of its past (and present) successes. It is in this sense that 
modern cosmology should be regarded as successful, and this will remain in the future, even 
if it might be regarded as being "wrong"... 

^ The choice of the coordinates system lead to a form of the metric for which the coefficients 
arc constant. 



proposed that this effect was th e ori gin of known redshift and provided the first esti- 
mation of the Hubble constant l|l41l) . In 1929, Robertson pubhshed the now so called 
Robertson- Walker metric. In the same paper he mentioned that distant sources appear 
to have a frequency shift. 



2 Basics of Friedmann-Lemaitre Models 

The fundamental idea of the geometrical theory of gravity starts from the fact that 
we can assign four coordinates to any event observed in our vicinity, for instance in 
Cartesian coordinates : {x,y,z,t). Locally, space appears flat to be. However this does 
not determine of the geometry of space at larger scales : local observations put us in 
the same situation that led people to think the earth was flat : the fact that we can 
describe our vicinity by a flat map does not determine the actual geometry on larger 
scales. Let us take the line element of a homogeneous 3D space which can be showqj 
to be : 

dl"^ = r'^ide'^ + sin^ 9d(l>^) + ^^— — 2 (2) 
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where k is —1,0, 1 according to whether space is hyperbolic, flat or spherical. -R is a 
characteristic size (in the spherical case, that is the radius of the 3D-sphere embedded 
in a 4D space). 

We then add the time as the fourth coordinate to build the equivalent of the 
Minkowski space-time element of special relativity and get the Robertson- Walker (RW) 
line element after the change of variables -Sr — >■ r: 

ds'^ = -C^dt^ + R(tf[r^(de'^ + sin^ 9d<l>^) + '^^ J (3) 

1 — fcr^ 



2.1 Topology 

The above line element describes the local shape of space : the curvature is only a local 
property of space, but does not tell us about the global shape of space. For instance, the 
Euclidean plane is an infinite flat surface while the surface of a cylinder is a 2D-space 
which is flat everywhere but is finite in one direction. GR will in principle allow us to 
derive the local geometry of space and its dynamics, but does not specify of the global 
topology of space. Only direct observations would allow to test what the topology 
actually is. Of course this will not be possible on scales much larger than what can be 
observed (the horizon). We can therefore hope to prove that the Universe is finite, if it 
is small enough, but we could not know whether we are in a finite Universe of which 
the scale is larger than the horizon, or whether we are in an infinite Universe. The 
interest in the topic of the cosm ic topology, with possible observational signature, has 
been recently revived l|l37l : 11521 ) 



•^ It is an instructive exercise to start from an Euclidean 4D space x,y,z,u and derive the 
line elements dl^ on the 3D sphere (x^ -(- j/^ -f z^ -|- m^ = R^) in internal spherical coordinates 

(r = ^a;2 + y'^ + z^, 0, <f)). 



2.2 Dynamics 

The function R{t) which appears in the RW hne element, is totally independent of any 
further geometrical consideration. It can be specified only within a theory of gravity. 
The basic equation of GR relates the geometrical tensor Gij to the energy-momentum 
tensor Ty 

G^j^R,j-^g,jR = 8-KGT,j (4) 

where gij is the metric tensor, Rij is the Ricci tensor, R the Ricci scalar. For a perfect 
fluid, there exists a coordinates system, called the comoving coordinates, in which the 
matter is at rest and the tensor Tij is diagonal with Tgo = p and Tn = T22 ~ T33 = p, 
p being the density and p the pressure. A fundamental aspect of GR is that the source 
of gravity includes explicitly a term coming from the pressure : p + 3p/c . Finally, there 
is an analog of the Gauss theorem, that is the Birkhoff's theorem (|26n I: if the matter 
distribution is spherical then the evolution of the radius of a given shell of matter 
depends only on its content. 

From the above rules, we can easily derive the equation for R{t). Let us consider 
a spherical region of radius a in a homogeneous distribution of matter. The equivalent 
Newtonian acceleration is : 

with the acceleration being generated by the "mass" M{a) of the above spherical 

region : 

GMja) 4 2^ .„s 

5 = -^ — ^ --^-^Gip + Sp/c )a (6) 

The density term includes the effect of kinetic energy {E = mc !). Writing total energy 
(Et) conservation inside the volume of the sphere from elementary thermodynamics 
gives : 

d{Et) = d{pVc^) = -pdV (7) 

leading to : 

P = -3(J + p)^ (8) 

From these two equations, the pressure can be eliminated, and, after having multiply 
both terms by d, the differential equation can be easily integrated. This leads to the 
following equation : 

'<i\2 SnGp Kc^ 



3 a2(t) ^^) 

The last term corresponds to the constant of integration. Its value cannot be specified, 
depending on the initial conditions. The form of the above equation is independent of 
the radius a of the sphere and the solution a{t) should be proportional to the quantity 
R{t). R{t) should also be solution of an equation of the same form, the constant K, 
which depends on the radius aq, being related to the constant k which is involved in 
the Robertson- Walker metric element, something which can be established only within 
GR: 

R\^ _ 8-kGp kc^ 



Rj 3 R^{t) 



* Apparently, this theorem should be named Birk hoff- Jebsen, as it has been published two 
years earlier by an Norwegian physicist, J.T. Jebsenl ll21h . 



This relation expresses the hnk within the framework of GR between the geometry 
and the material content of the Universe. In order to specify completely the function 
R{t), one needs an equation of state for the content of the Universe. The two cases 
often seen in cosmology are the dust case (p — 0) and the radiation dominated regime 

(p=ipc2). 



2.3 Vacuum and the Cosmological Constant 

Vacuum is a particular medium, and one should wonder what is the equation of state 
of this medium. Naively, one would think that the equation will be pv = and pv = 0. 
However, let us try to derive the vacuum equation of state from first principles. As in 
classical thermodynamics let's assume that we have a piston with vacuum in it. We 
also assume that simple vacuum (p = 0, p = 0) is present outside. 

The energy inside the piston is _B = pvC V. If the volume changes by a small 
amount the net energy change is: 

dE = dip^Vc^) = pvC^dV (11) 

this change is equal to the work of the pressure: 

dE = -pvdV (12) 

so the equation of state is: 

Pv = -pvc (13) 

As one can see, the conditions p > and p > ensure that the simple solution is the 
only one. However, there is nothing which imposes these conditions for the vacuum, 
and we can therefore decide to keep such a possible term. This can be directly translate 
in the equations governing R{t) by introducing the following constant : 

A = SttGpv (14) 

Such a term is called the cosmological constant and has been historically introduced by 
Einstein as a modification of his original theory. It appears as an additional term in the 
left hand side of eq. |4] We have recovered the two usual Einstein- Friedmann-Lemaitre 
(EFL) equations: 

rV _ StvGp kc^ A 
and 



_R _ 47r_G (^ , 3p\ , A 
R ^ 3 



(.+5)+ 1 («) 



There are several ways to write the EFL equations. For instance we used energy con- 
servation, eq. ISl which can be kept instead of one of the above eqations. It can also be 
useful to use: 

H = -4.g(p+J)+^ (17) 

Only two independent equations were obtaineded, while three unknown quantities are 
to be determined {R{t) , p{t) , p{t)) , therefore we need a further constraint which is 
provided by the equation of state F{p,p) = 0. 



It is quite usual to write the fundamental cosmological quantities as non-dimensional 
quantities which depend on redshift. The following notations are very common : 
H — %, the Hubble parameter, 
Oj^l — fl = -TTrr the density parameter, 

_ _ii^ the deceleration parameter, 
Qv ~ fi\ — \ = ^TWi, the (reduced) cosmological constant, 

, 2 

^c = —jTrnj ~ ~<^j the curvature parameter. 

Quantities are labeled with a when they are referred to their present value. For 
instance the present day value of the density parameter f2 is Oq. With the above 
notations, the first EFL equation 1 151 then reads : 

^M + nc + n^ = i (18) 

or: 

a = n + \-l (19) 

so that the "radius of the Universe" can be written : 



while the Hubble constant evolution is given by: 



i/2 = H^(J7m(1 + zf + Oc{l + zf + n^) (21) 



2.4 Important quantities needed for observations 

In this section we only need to work in the framework of a geometrical theory of space- 
time, in which the trajectories of light rays are assumed to be the null geodesies. Let 
us have a comoving spherical coordinate system (r, 9, (p, t) the observer being at the 
origin of the spatial coordinates (r = 0, ^ = 0, = 0, io)i let assume that the observed 
source is emitting light at the coordinates (rg,6l = 0, (/) = 0,ts)i a^nd let r{t) be the 
trajectory of the emitted photons. As this trajectory is a null geodesic, we have: 

'^'dt'^R\t)^!^ = Q (22) 



so the variables can be separated and the integration over r is analytical: 



Its 
with: 



( sin(rs) if fc= +1 
Sfc(r5) = <^ rs if fc=0 (24) 

I, sinh(rg) if k — -1 

When the distance is small with respect to Rq we just have S7 ir) ~ r. 



2.5 The Redshift 

In order to derive the observed frequency uq of the hght from a source emitted at the 
frequency v, we consider the trajectory of a second light ray emitted at the time is + ■^• 
As the source is comoving its coordinate is unchanged and we have: 

1 /*« cdt ['"+'/'"'> cdt 



(26) 



which implies: 

1^0 ^ Xs ^ Rs ^ 1 

V Aq Rq 1 + z 
where z is the redshift. This is the standard formula for the cosmological shift of 
the frequencies. This result shows that the redshift 2 is a natural consequence of the 
expansion. 

2.6 The proper distance 

In GR, space changes with time, and there is no proper time, so that the "intuitive" 
notion of distance between two points is not a well defined quantity. Therefore the 
various methods to measure the distance between an observer and a given source give 
different answers. The proper distance - between the source and the observer - can be 
seen as a distance measured by a set of rulers at time t. The distance element is given 
by: 

dl^ ^ ds^ = Ritf '^ ^ (27) 

so that the proper distance is : 

Dp ^ R{t)S^\rs) (28) 

The fact that this distance changes with time is the direct consequence of the expansion 
of the Universe. We can now examine how this length changes with time : 

Dp = RS-\r) (29) 

so that the source is actually receding from the observer with a speed: 

V^Dp^ ^Dp = HDp (30) 

ri 

The fact that this speed could be larger than the speed of light should not be considered 
as a problem: this speed can be measured but cannot transport information faster than 
light. When the distance is small, the Doppler frequency shift is : 

'-^ = ist = H^ = ^ (31) 

u R c c 

so that the shift is the one corresponding to the Doppler shift associated with the 
above velocity. For large distances, the total shift results from the product of small 
Doppler shifts and the redshift is therefore purely kinematic. The ph ysic al natu re of 
the expansion has been recently the subject of interesting discussions (|5ll: 111: 11761 : ISTh . 



Comoving distances It is sometimes useful to refer to comoving distanceqj. The co- 
moving distance D'^{z) associated to the distance D{z) is : 

D'{z) = ^D{z) = {l + z)D{z) (32) 

In the case of the proper distance, this becomes: 






*° cdt r dz 



R{t)/RQ 7o H(z 



(33) 



2.7 The angular distance 

Let us suppose that we observe a ruler orthogonal to the line of sight. The extremities of 
the ruler have the coordinates (r, 0, 0, tg) and (r, 6, 0, ig). The proper length I between 
the extremities is: 

l^ = ds^ = Ritsfr'^e'^ (34) 

which provides the relation between the angle 9 and the length I and thereby the 
angular distance defined by: 

7?a„g = !^= Rits)r (35) 



2.8 The luminosity distance 

Let us assume that we observe a source with an absolute luminosity L through a 
telescope with a diameter d and let us choose a coordinates system which is centered 
on the source. Let 9 be the angle between two rays reaching two points diametrically 
opposite on the telescope. We have d = R{tf))r9. The energy emitted by the source 
that reaches the telescope is : 

When observed, the energy of photons has been shifted by 1/(1 + z) but also the 
frequency at which they arrive is reduced by the same factor. Therefore the flux (energy 
per unit time and unit surface) one gets is: 

f- ^ 1 - ^ - ^ (37) 

•' 7r/2/4 (1 + 2)2 47ri?(io)2r2(l + 2)2 47rD2^^ ^ ' 

This relation provides the luminosity distance: 

Aum = R{to)r{l + z) = i?(ig)r(l + zf = Dang(l + zf (38) 



2.9 Distance along the line of sight 

We consider here the length along the path of a photon trajectory. The length element 
is 

dl = cdt^c^ = --^-^ (39) 

R H{z)l + z 

This relation is useful to write the volume element. 



This could also bo confusing! 



2.10 The age of the universe 

The general expression of time interval is 



7i d-R 1 dz ,^„, 

R H{z)l-\'Z 



3 Some solutions of the EFL equations : relativistic cosmological models 

The various possible theories of gravitation provide different functions R{t) and through 
the above tests may in principle be distinguished by observations. However, it is easy 
to check that the difference only occurs at high redshift. In practice these tests may not 
be discriminant because the observations of distant objects are difficult and because 
the universe at high redshift is younger, so any object at high redshift is likely to be 
different from to-day because of evolution. In other words, general geometrical tests 
rely on the assumption that the typical evolution time scale of the objects undr study 
is much larger than the age of the universe. 

One must also underline that as already mentioned the EFL equations can be solved 
only once the equation of state is specified, i.e. a relation between p and p is adopted. 
It is now quite common to specify this relation by the w parameter: 

p — wpc (41) 

If w is constant, then eq.|8] allows to find the evolution of the field density: 

p{z) = po{l + zf'-'+-^ (42) 

while for a general function w(z) one has: 

piz) = poe^pfJ\{l + w{z))^\ (43) 

There are three particular regimes : the matter dominated one, the radiation dom- 
inated one and the vacuum dominated one. 

In the matter dominated case one has p = 0, so ui = and the mass (per comoving 
volume) is conserved : 

pa^ = p{l + zy^ ^ cste (44) 

while in the pressure dominated case p = -^pc , so w = -^ 

pa = p{l + z) — cste (45) 

Finally in the vacuum dominated case p = —pc , so w = —1 and 

p = cste (46) 

In some models lu < — 1, and therefore the density of the universe increases when it 
expands. 
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3.1 Case (Aq = 0,p = 0) 

When the cosmological constant is zero, there are three types of solutions: 

a) when the density is above the critical density: 

P>Pc = |^ = 2.10-^Vg/cm3 (47) 

the spatial solution is the spherical space. The function R{t) grows from zero to a 
maximum value then a collapse phase follows to zero. 

b) when the density is equal to the critical density, the solution is named the 
Einstein-de Sitter universe and R{t) is simple : 

R{t) = Ro [^Hotj = Ro [-) (48) 

with to = ^H-^ = l/(67rGpc)^/^ 

c) when the density is below the critical density, the function R{t) grows from zero 
to infinity. It is easy to check from the EFL equations that the function R{t) behaves 
like t when R is large. 

The behavior of R{t) can be found when t — )■ independently of the model: R{t) oc 
i2/3. 

Finally, the relation between the comoving coordinate r and the redshift can be 
expressed: 

_ c 2 i?o(l + z) + 2-2J7o -(2-r?o)Vl + J7oz 
°'^ ' Hon^O 1 + z 

This is known as t he M attig relation. Others useful quantities can be found in Weinberg 
recent text book (|233l ). 



3.2 Case (Aq >0,p = 0) 

There are many possibilities when a cosmological constant is allowed. To specify a 
cosmological model, it is customary to specify two " observables" : i?o ^-nd f2\. For 
instance, the cosmological view of the Friedmann-Lemaitre models is summarized in 
figure \T\ 

The look back time, i.e. the time since the epoch corresponding to the redshift z 
con be obtained from 1401 

if du 



^(2) = -FT / , , , = (49) 

^0 Ji uy i?o^ ~ aQU^ + Aq 

where u = 1 -[- z. When A < it can be shown from the expression of R that 

to = r((») < -^ (50) 
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Fig. 1 Evolution of the scale-factor R{t) with time, in FL models according to the values 
of the parameters i7jv/ and f2yi . One notes that above the critical line universes has no initial 
singularity; on the right side of the second branch of the critical line are the re-collapsing 
models. 



3.3 Radiation dominated case 



As we have seen the density associated with radiation evolves according to p^a = cste. 
It is clear that this radiation term will be dominant over the matter term for small 
a, that is at the "beginning". In this regime, if we neglect other terms in the EFL 
equation, we have: 



R^ 



V 3 y R 



(51) 



so that Roct 



1/2 



3.4 Basics of structure formation 



Structure formation in an expanding Univ erse through gravtational instability can be 
traced back to the early works of Lemaitre (|l42r ). who established the spherical solution 
and discussed its relevance to cluster formation. In modern cosmology, linear and non 
linear regime {5pp 3> 1) are used to obtain useful constraints. 



The gravitational instability picture The linear regime of small perturbations (both in 
size and amplitude) in the matter field (p{x) = p (1 4- &{x))) in the absence of pressure 
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can be derived from newtonian equations (|17' 

S + 2HS = ATvGpmS = |fim-ff^<5 (52) 

Because there is no spatial derivates in the equation, the evolution with time is in- 
dependent of the scale (this does not hold any more when the pressure is taken into 
account). There are generally two distinct solutions to this equation, one is a decaying 
mode and the second one is a growing mode which is of relevance for cosmology. The 
growing mode is usually written as D{t) with a normalization such that the present 
day value is one: Dito) = 1. The evolution of a perturbation can then be written: 

m = D{t)So (53) 

where Sq is the linear amplitude that the fluctuation would have today. The qualitative 
behavior of solutions of eq. [S2]is easy to understand: if Om ^ 1 the right hand side 
is zero and the growing mode is frozen, while if Om ~ 1 there is a growing mode 
with D(t) oc a(t). The form of eq. 1521 implies that once the expansion rate history 
{H{z) and the present day density of the universe po ^re known, the solution D{t) is 
entirely determined (the matter density parameter can be inferred from its definition 

fim(z) = ^ = fio „ ). Therefore having a measure of the growing rate of 

3H-' H[z)-' 

perturbations is a direct test of the relevance of Newtonian gravity on the scale under 
consideration. Theories beyond GR might lead to growing rates that depend on the 
scale and departs from their standard values. 

Non baryonic dark matter and Cold Dark Matter scenario In the early universe the 
growth rate of fluctuations is more complex : pressure has to be taken into account, 
and there are at least two matter contents to take into account: baryons and the 
non-baryonic matter. These ingredients will introduce scales in the dynamics of fluc- 
tuations. Ultimately these scales will leave visible imprints in observables quantities : 
the fluctuations of the cosmic background radiation and the power spectrum of the 
matter distribution (or equivalently in the correlation function). The first scale comes 
from the difference in the dynamics for fluctuations larger than the horizon compared 
to the ones that are smaller. Indeed, when photons (and other relativistic particles) 
dominated the density of the universe the right hand side of eq. [52] is negligible, and 
therefore the growing mode is frozen. However, on scales larger than the horizon the 
dynamics should be evaluated within a full GR treatment. Intuitively, on these large 
scales the radiation acts as a source of gravity which will allow the fluctuations density 
to grow. Therefore a scale appears corresponding to the scale of the horizon. When the 
universe enters the matter dominated area, fluctuations on all scales will grow provided 
that the pressure term is negligible, i.e. for the non-baryonic matter. The scale which 
is imprinted is therefore the size of the horizon at the time of equality between matter 
and radiation and is characteristic of Cold Dark Matter (CDM). However, the baryons 
are still tightly coupled to the photons and thereby their dynamics is under the control 
of the pressure term to be added in the rigt hand side of eq. 1521 This coupling stops 
when the matter becomes neutral enough (z ~ 1500). A second scale is therefore im- 
printed in the baryonic component corresponding to the scale of the horizon of sound's 
velocity at the epoch of recombination. Although the baryonic matter is subdominant, 
the distribution of dark matter will nevertheless reflect this scale. This feature appears 
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as oscillations in t he power spectrum, but more distinctly as a peak in the correla- 
tion function (|160l ). In practice, the dark mater distribution is described by its linear 
power spectrum Pdm (k) at present epoch corresponding to the initial power spectrum 
Pi{k) modified by the evolution. In the linear regime the various wave numbers are 
independent so: 

PoMik) = Pr{k)TDM{k) (54) 

The nature of the initial fluctuations, the nature of dark matter (hot, warm, ...) are 
other ingredients that determine the shape of the transfer function. In practice, the 
computation of the transfer function is to be done numericall y. In the cold dark matter 
case (CDM) analytical functional forms have been provided (11781 : 1 14 : ITJ ) including the 
detailed baryonic features (J72n . 



4 The Quest for cosmological parameters 

4.1 Hubble diagram 

Looking for a possible relation between apparent velocities and distances was a scientific 
question that few astronomers addressed before Hubble's discovery. The determination 
of Hubble's constant have been the subject of considerable effort over the next eighty 
years following his discovery. One of the major scientific question that the Hubble 
Space Telescope was aimed to address was precisely the determination of the distance 
of close galaxies and thereby to offer an accurate measurement of the Hubble constant. 
The final value obtained from the HST key project to measure the Hubble constant 
concluded Hq — 72±8 km/s/Mpc |SJj), also the debate is not entirely clo sed a s Sandage 
and Tammann still concluded to a lower value: Hq — 62.3±4 km/s/Mpc l|214f) although 
these values are marginally consistent given the quoted - systematic - uncertaintieo 

4.2 Going further and other geometrical tests 

The first cosmological tests which have been proposed were based on the relation 
between some observable quantity to a corresponding intrinsic property of the source 
(like the relation between the apparent luminosity and the intrinsic luminosity). These 
tests are essentially geometrical in nature, they involve the coordinate r given by eq. 
l25l or the relation between the coordinate r and the redshift z: 

the dependence on cosmological parameters coming from the difference in the expansion 
rate according to eq. [TJl The Hubble diagram which has been extended to high redshift 
thanks to type la supernovae is the most popularized example of these geometrical tests. 
It is also rather intuitive that the Hubble diagram which expressed the speed of the 
expansion will provide direct information on the rate of acceleration/deceleration when 
extended to objects distant enough to be seen at appreciably early epoch of the history 
of the universe. A reliable extension of the Hubble diagram to high redshift has been 
made possible thanks to the use of type la Supernovae. SNIa could have a maximum 



The Hubble constant is traditionally noticed: Hq = hlOO kni/s/Mpc. 
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luminosity (M ~ —19.5) comparable to that of an entire galaxy. Furthermore there is 
a relation between the decline rate and the intrinsic luminosity making them suitable 
for distance measurements at cosmological scale. Because SNIa are rare, large sky area 
have to be surveyed on a regular basis to collect samples of SNIa. At the end of last 
century, two groups have independently investigated the distant SNIa Hubble diagram 
and concluded that supernovae at redshift ~ 0.5 were dimmer by ~ 0.2 mag compared 
to what was expected in a unaccelerated universe. This was interpreted as an evidence 
for an accelerated expansiorjj ( . 196: : 186 : 101: 85) . 



4.3 The age of the Universe 

Formally, we have no direct information on the actual age of the universe to- But there 
are several astrophysical objects for which an age can be derived. The most common 
one is probably the age of globular clusters. There are different limitations which 
make determinations of absolute ages problematic. Distance estimations is for instance 
one limitation. Hipparcos results have helped, but the uncertainty remains significant. 
Recent estimate for the age of globular clusters was given to be 12, 6T3'4(95%) Gyr 
(|l34i '). an age consistent with that of the star CS 31082-001 estimated to be 14. ±2.4 Gyr, 
based on the decay of U-238 ([109|). It has been noticed since the beginning of cosmology 
that age estimations were high given the preferred value of the Hubbl e con stant. This 
has been noticed by Lemaitre when debating with Einstein, Peebles (|l81h and many 
others since, providing thereby an argument in favor of the cosmological constant which 
holds until the late 90's (ll35J;[7j). However, this argument is relatively weak because 
of the difficulty to achieve reliable astrophysical estimates of -ffo- With modern best 
values, the product HqIq can be evaluated numerically: 

Hoto « 0.102/ltGyr (56) 

with h — 0.72 one gets HqIq ~ 0.93 while only values greater than 1 would request 
an accelerated universe. Although, it is important to have a consistency check, im- 
provement by one order of magnitude in precisoin will be necessary in order that age 
and Hubble constant values could be combined to provide constraints on cosmological 
parameters competitive with current constraints. 



4.4 The Cosmic Microwave Radiation 

Spectrum and uniformity The measurement of the CMB spectrum to check for the 
black-body shape of this spectrum was one of the most important cosmological test 



^ There is an intense debate on the "discovery" of acceleration from SNIa Hubble diagram 
in order to specify who did what and who said what. The hist ory o f the research program 
developed within the SCP is summarized by Gerson Goldhaber ll lOOl l . An early popular sci- 
entific report on SCP talk at AAS meeting January 1998 meeting has been published by 
Jame s Glanz (Science, 279, 651), while some view from the High Z team by Robert Kirshner 
il26l ) is available at http://www.cfa.h5irvard.edu/~rklrshii6r/whowhatwhen/Thoughts.htni 
It is also interesting to mention that claim for evidence of an accelerated expan- 
sion from a totally different technique has been presented at the same AAS meeting: 
http://www.bk.psu.edu/faculty/daly/PU98.pdf These data taken at face value were point- 
ing to a low density universe, but only the flatness of the universe as already evidenced by 
CMB fluctuations would have allowed to conclude on the actual acceleration of the expansion. 
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Fig. 2 The amplitude of angular fluctuations of the CMB is expressed through their angular 
power spectrum. Data are WMAP, Boomerang, ACBAR (193!). A simple minimal six parame- 
ters model including a cosmological constant provides an excellent fit to the data. This is one 
of the most important successes of modern cosmology. 



awaited for since the discovery of the backgroun d radiation. Two high quaUty mea- 
surements have been pubhshed in 1990 ( 1051 : 11591 ). while the FIRAS final results were 
remarkably accurate (86): any departure from a blackbody should represent only a tiny 
perturbation of the total energy: 5E/E < 10~ . 

An interesting application of cosmic background radiation is to provide a test of the 
reality of the expansion: if we are able to measure the temperature of the background 
at higher redshift it should scale according to: 



r(z) = To(i + 2) 



(57) 



It is actually possible to measure the temperature of the background radiation through 
the observations of the ratio of molecular lines: the ratio of the population on two 
levels for which the difference in energy is only a few Kelvin is sensitive to the cosmic 
radiation field and therefore provides a way to actually measure the temperature of 
the background. Such lines can be detected in the optical domain. Actually the first 
detection of the CBR was obtained by this method with CN lines (162 : 241 ) . It has also 
been successfully applied to distant QSO's to measure the variation o fthe temperature 
of th e background with redshift and the result are consistent with the expanding picture 
(|21l|;[213). 

These measurements provide a fundamental test of the Big Bang picture but do 
not provide a source of constraint on the cosmological parameters. 



Fluctuations Since the discovery of the CMB fluctuations by COBE (l210l) the idea 
that early universe physics has left imprints revealed by these fluctuations has gained 
an enormous attention. In this respect, DMR results have played a fundamental role 
in modern cosmology comparable to the discovery of the expansion of the universe or 



16 

the discovery of the microwave background by Penzias and Wilson, and indeed this has 
motivated the dehvering of the Nobel prize to G. Smoot and J. Mather for this discov- 
ery. One of the fundamental reasons for this is that fluctuations on scales larger than 
one degree in the microwave background radiation correspond to scales greater than 
the horizon at last scatering epoch and cannot therefore b een a ltered by any physical 
process and should therefore reflect primordial fluctuations (|23in . This also means that 
the very existence of these fluctuations could be explained only from yet undiscovered 
physics, probably relevant to the very early universe |10Q), for which the expansion law 
has been modified compared to the standa rd p icture. The DMR results were providing 
some constraints on cosmological models (242) but it has been realized that the mea- 
sure of fluctuations on smaller scales will provide much stringent information. Early 
detection of fluctuations on degree scales allowed to set interesting constraints and 
provided the first evidence for a fiat geometry of space (|l47l : ll48l ). If estimations of low 
matter density were to be regarded as robust, this was inevitabily leading to a non-zero 
cosmologcial constant. Even before the availability of the WMAP data, considerable 
progresses have been achieved on the measurement of fluctuations on all angular scales. 
Archeops (1201 ) and Boomerang (|22h . as well as many other small scale measurements, 
already provided measurements allowing tight constraints on cosmological parameters 
(|2ll ). Although the observed fluctuations were consistent with a A dominated universe, 
a cosmological constant was not explicitly requested by the CMB data alone. Indeed 
even the WMAP data were consistent with a vanishing cosmological constant, provided 
the Hubble constant was left as an entirely free parameter. A positive detection of a cos- 
mological constant could be obtained only by using some additional data in conjunction 
with CMB [21 like the measurement of the Hubble constant. A further restriction came 
from the fact that the constraint on cosmological parameters were obtained within the 
standard CDM picture, and that many ingredients were specified without being nec- 
essarily confirmed by observations : for instance initial fluctuations are sup posed to be 
adiabatic and to follow some power law. Therefore the "concordance" (|l73h cosmology 
was an appropriate terminology: the model was consistent with most existing data, but 
the introduction of a cosmological constant was not requested by any single data, and 
it was far from being clear whether relaxing some of the input hypotheses would not 
allow for solutions without the introduction of a cosmological constant. 



4.5 Clusters of galaxies 

Most talks on clusters of galaxies start with stating the fact that they are the largest 
virialized structures in the universe. Another reason for the strong interest they repre- 
sent is due to the fact that clusters are the only astronomical large objects for which we 
have such a wide range of information: their total mass is in principle accessible from 
x-ray spectroscopy, velocity dispersion, weak lensing signal, their gas content can be in- 
vestigated through x-ray observations but also through their imprint on the microwave 
background, the thermal Sunyaev-Zel'dovich effect. Obviously their star content can 
be evaluated, but also the metal content of their gas. Clusters are therefore the astro- 
physical objects relevant to cosmology on which the most comprehensive set of data 
can be obtained. 
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4-5.1 Baryon fraction 

X-ray data allow to determine at the same time the mass of the gas content of clusters, 
through imagery, as well as the gas temperature., through spectroscopy. The temper- 
ature can be used to estimate their total mass. It is therefore straightforward then 
to infer their baryon fraction (most of the baryons being in the gas, the star even if 
uncertain can be taken into account without making much difference). Because, the 
gravitational collapse is expected to be identical for dark matter and for baryons, at 
least before the first shocks occur, the baryon fraction should reflect the primordial 
baryon fraction : 

where T is expected to be close to one. Numerical simulations shows indeed that T 
should lie in the range 0.85 — 1.. i?^ can be known from primordial nucleosynthesis 
results, and therefore constraints can be obtained on SIm l|235f l. Although this test 
is attractive in its principle, its application, especially in the context of precision cos- 
mology, is somewhat problematic. There are several reasons for that: for instance, in 
numerical simulations, as well as in the observations, T seems to depend on the radius 
and is therefore a function T{r) on radius. The precise value of T{r) is likely to depend 
on the physics of the non-gravitational heating and cooling processes which are nec- 
essary to reproduce observed properties of clusters in numerical simulations. However, 
these possible limitations still allow to infer an upper limit on i7jv/ from eq. 1581 A final 
limitation comes from the clumping of the gas which would lead to an underestimation 
of Om- On the observational side, it is unclear to which accuracy cluster masses are 
known : agreement between lensing masses and x-ray mass is under debate. Estima- 
tions of total mass and baryon mass in the outer parts of clusters, which are expected 
to provide the most reliable estimation of the primordial baryon fraction, are subject to 
uncertainties: the x-ray emission of the gas at the virial radius is low and only for few 
clusters do we have a firm detection. At this outer part, the gas is eventually clumped 
which would bias the gas mass estimates. 

4-5.2 Geometrical tests 

Another strategy is to use clusters as standard candles in some sense and to constrain 
the cosmological parameters from the observed properties of distant clus ters. For in- 
stance it has been proposed to use cluster x-ray radius as a standard ruler (|l67h . Addi- 
tional information may allow direct distance measurements. This is what the Sunyaev- 
Zel'dovich (SZ) signal provides: because the x-ray emissivity goes as the square of the 
gas density, while the SZ signal just goes as the density the distance can be directly 
obtained. Not only this method can be used to estimate the Hubble constant, it is also 
possible to extend the Hubble diagram to distant clusters and thereby to get the same 
information provided by the SNIa diagram. Another variant of these tests has been 
proposed from the baryon fraction of clusters against redshift: because the gas fraction 
inferred from observations depends on the angular distance to the cluster, requesting 
that the gas fraction does not evolve with redshift will result on constraints on the 
cosmological parameters (20S). This test has deserved recent attention (jj), but uncer- 
tainties mentioned about the baryon fraction are l i kely to alter this test either and it 
is difficult to draw robust conclusions from it (|202l :l83l:l79l). 
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4-5.3 The Use of Cluster Abundance 

The theoretical mass function A theoretical expression for the abundance of virialized 
structures in hierarchical picture has been first proposed by Press and Schechter in 
1974 y_89j}. From scaling arguments, a general expression for the mass function can be 
derived©: 

Snl{z) being the threshold for non-linear collapse of spherical density perturbations 
(in an EdS universe 5ml{z) = 1-68), az{M) the amplitude of mass fluctuations on 
the mass scale M (ctz(M) = ao{M)D{z), D{z) being the normalized growing rate of 
linear fluctuations) and vml the normalized threshold {vml = Snl/'^{M) ) for non- 
linear collapse, J T{v)dv is the fraction of space covered by spheres of mass > M 
satisfying the non-linear criteria. Such an expression can be rigorously justified only for 
a given power law spectrum in an Einstein-de Sitter universe. However the numerical 
investigations of the mass function have shown that it does follow the above scaling to 
a very high accuracy (|69i : i207i : illS) . Press and Schechter used for the function J-"; 

2 . y"" 



.FH = ^/-exp(- — ) (60) 

TV I 



while the expression proposed by Steh and Tormen (|207l : 12081 ). possibly extending the 
non-linear condition to take into account ellipsoidal collapse, was found to produce a 
more accurate agreement with numerical simulations : 

^{v) = a/ — Cexp(-0.5A//2)(1. + {l./{Auf)'^) (61) 

y TT 

with A = 0.707 C = 0.3222 and Q = 0.3. Jenkins et al. (flli l did provide another fitting 
formula in which the amplitude of (5jvl is set constant (independently of the cosmology) 
which was found to provide an accurate fit to within 20% in relative amplitudqj. This 
accuracy is sufficient for present-day available data to be compared with. More recently 
larger numerical simulations were used (|23d : 12171 : 1591 ) to investigate the mass function 
to a precision higher than 20%, finding some departures from strict universality and 
provide more accurate fitting formula, but cosmological data are not yet at a precision 
level for which such departures could be observed. 

Local abundance The observed local abundance of clusters of a given mass M can be 
used to infer directly the average amplitude of matter fluctuations, a{M), through the 
expression of the mass function. Such a derived value depends weakly on Cosmology 
(JSOJ) as Cosmology enters only through iStvl, which as we have mentioned is almost 
independent on Cosmology. However, the traditional way to express the amplitude of 
matter fluctuations is through ag and therefore the dependence of this quantity appears 
on Cosmology and on the initial spectrum. There is a source of uncertainty due to the 
fact that the mass is not directly an observable quantity. Therefore one has to use a 



* Another complication happens from the definition of "an object with a mass M" . Two 
algorithms are commonly used, the friend of friend one (FOF) and the spherical overdensity 
one (SO). In addition the non linear contrast density A which is used could be set to an 
arbitrary fixed value or to depend on Cosmology according to the virial value Av evaluated 
within the spherical collapse model. 
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relation between mass and some observable quantity. At the same time, the selection 
function of clusters within a given sample has to be known in order to infer the actual 
number per unit volume. X-ray clusters are a priori the easiest choice, as the selection 
function from x-ray surveys is easy to model and the x-ray temperature is a priori a 
reliable and physically motivated proxy of the total mass. There are various arguments 
from which one expects that Tx oc GM jr. This leads to a scaling relation between 
mass and temperature: 

T = Ata/M^^/^(J7m(1 + -4)/179)^/^/i^/^(1 + 2) keV (62) 

an expression in which the contrast density A is the relative density compared to the 
background density of the universe. An other density contrast is also used which refers 
to the contrast density Ac relative to the critical density of the universe pc(^z). The 
above scaling reads then: 

T^ATMMl^^{nM(X^Ac)l\l'dil^{E{z)hfl^ keV (63) 

where the redshift dependent Hubble constant is written Ii{z) = 100/i-B(z)km/s/Mpc. 
There has been much debate on which val ue of the normalization Aj^m is to be used 
and this translates in an uncertainty in erg (|l88J). The various values published for erg, 
ranging from 0.6 to 1.05 for a ylCDM model with fl = 0.3 l|32r ) are consistent with 
each other when the various values of the normalization constant ^ta/ a^re taken into 
account. However, because of this uncertainty t he a ctua l val ue of ag inferred from 
cluster abundance is still subject to some debate (|l95l : l8ll : ll97l ). 



Abundance evolution Peebles et al. (|180f ) and Evrard (|80f ) discussed the constraints 
that the existence of high redshift clusters were setting on Cold Dark Matter models, 
at a time when Einstein-de Sitter Cosmology was fashionable. These works followed 
some earlier ones on the use of cluster abundance in cosmology ([187|). The evolution 
of the abundance of clusters of a given mass is a sensitive functi on of the growing 
rate of fluctuations therefore offering a powerful cosmological test ([14J; [l7J). Rather 
obviously, this test needs that local abundance to be adequately normalized. This test 
is primarily sensitive to the cosmological density of the Universe, very weakly depend- 
ing on other quantities (|30f). Attempts to apply directly the test of the evolution of 
the temperature distribution function of clusters have been performed but still from a 
very limited number of clusters (typically 10 at redshift 0.35) and lead to somewhat 
apparently conflicting results (|108i : i226 : i31il . 

On the other hand, redshit number counts (without temperature information avail- 
able) allow one to use samples comprising many more clusters. Indeed using simulta- 
neously different existing surveys one can use information provided by more than 300 
clusters with z > 0.3 (not necessarily independent). In order to model clusters number 
counts, for which temperatures are not known, it is necessary to have a go od knowledge 
of the L — T relation over the redshift range which is investigated (|l75|). Early inves- 
tigations indicated the absence of evolution of t he L — T relation or a slight positive 
evolution, indicative of a high density universe (|20ll : 1 194 : |38|). The L — T relation at 
high reds hift h as been more recently determined by the XMM-J7 project (|l6l : ll5H ) and 
Chandra (|228r ). Number counts can then be computed: 

7V(> U,z, 2Az) = n /;+^; ^{Lx > 47rDffx)dz 
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= " Itt ImTz) ^(^' ^)dMdV{z) (64) 

where T[z) is the temperature threshold at redshift z corresponding to the flux fx 
as given by the observations, being therefore independent of the cosmological model. 
This approach has been used to model all a vailab le ROSAT surveys and using the 
latter L — T measurement provided by XMM (12251 ). All existing x-ray cluster surveys 
systematically point toward high Om- A combination allowed a determination of Om 
with a 15% precision: 0.85 < Qm < 1.15(lcr) (depending somewhat on the calibration 
issue, this is part of the systematic). During this analysis numerous possible sources of 
systematics were investigated with great detail and are representing roughly an addi- 
tional 15% uncertainty. This means that global uncertainty is roughly 20%. This gives 
unambiguous evidence that the observed high redshift cluster abundance compared to 
local one is inconsistent with the one modelled within the concordance model. This is 
confirmed by the recently availabilit y of a large sample of high redshift clusters with 
measurements provided by Chandra (|229l 'l and a sample of high redshift massive clus- 
ters (|67l ). The temperature distributions drawn from these samples is shown in figure 
|3] : in agreement with previous results these temperature distributions are inconsistent 
with the concordance cosmology. A critical ingredient in the modeling of x-ray clus- 
ters is the mass temperature relation O A possible strong evolution of this relation, as 
proposed to reconcile redshift number counts with predictions (225), would allow the 
agreement of the concordance model with the observed temperature distributions as 
well. 



4.6 Large scale Structure of the Universe 

The properties of the galaxy distribution on large scales has been a major source of 
constraints for cosmology in the past. The most often quoted quantity is the correlation 
function of galaxies. Until the 1990's the correlation function was relatively well know 
only on scales which are non-linear or nearly so. Therefore the correlation function could 
be adjusted only by means of numerical simulations. Two major breakthroughs enabled 
important progresses in this area : the first one was the measurement of the correlation 
function of galaxy from the APM survey allowing to exclude the standard cold dark 
matter picture (153). Furthermore the data were consistent with the predictions of a low 
density CDM model with flf^j ~ 0.3 and it has been argued (|70|) that the introduction 
of a cosmological constant was a way to r econcile CDM with infiation predictions, 
in the spirit of an earlier proposition (|l79f ). The second one was the discovery of a 
simple analytical formalism allowing to construct the correlation function in the non- 
linear regime from the initial spectrum (jlOTI 'l. More recently, the most critical advance 
resulted from the availability of very large galaxy surveys, the 2Df redshift survey and 
the SDSS survey, allo wing to m easure the amplitude of galaxy fiuctuations on scales as 
large as 100h~ Mpc (|l83l : l215 : 73 : 184 ). This has provided a remarkable success to the 
ylCDM picture because the amplitude of the correlation function could be predicted 
for models that already match the CMB fiuctuations measured by WMAP: not only 
ylCDM model reproduces the shape of the correlation function, but the specific presence 
of a bump in the correlation function at scale of the order of WOh~ Mpc due to the 
detained dynamics of fiuctuations when the baryons are taken into account, the so 
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Fig. 3 The low and high redshift temperature distributions compared to predictions of an 
Einstein de Sitter Universe (left side, erg = 0.55, Atm = 6.35 in eq. , F = 0.135 is a shape 
parameter for the spectrum) and to predictions of the concordance model (right side, erg = 0.81, 
A-TM = 7.75, r = 0.168 ). Upper (blue) losanges are the local data. Lower (red) losanges are 
the temperature distribution function (z ~ 0.33)) derived from EMSS, (yellow) squares are 
high redshift (z ~ 0.5)) clusters from MACS and the 400 sq deg survey. 



called accoustic peak, corresponding to the "peak" in the C; of the CMB (Fig.|3|. More 
detailed use of pwoer spectrum measurements might be limited by our understanding 
of the exact rela tion between galaxies distribution and the underlying dark matter 
distribution looi). 



5 Possible origin of the apparent acceleration 

5.1 An Einstein-de Sitter Model ? 



ernovae Hubble diagram The first convincing evidence for acceleration is generally 
considered as coming from the SNIa Hubble diagram. However, use of geometrical tests 
based on the assumption of no-evolution of the parent population of the test is always 
possibly subject to produce biased values because of un-anticipated evolution. One 
possible way to cure this problem is to assume some evolution and see whether the 
data still provide evidence for the claim. For instance, an evolution term like: 



Anie 



(65) 



can not mimic the observed Hubble diagram. The limitation of this approach is due 
to the fact that we have no model for this evolution and therefore we are left with 
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Fig. 4 Large scale correlation function of the SDSS Luminous Red Galaxies ll73f) . Three 
ylCDM models (i7ji,fh^ = 0.12,0.13,0.14) are shown with a baryonic content evaluated from 
the CMB (fiijh^ = 0.024). A model without baryon corresponds to the bottom (magenta) line. 



a purely empirical approach. Other forms of evolution may therefore lead to different 
conclusions. Indeed, it has been suggested that the correcton term might be: 



AlUe oc At 



(66) 



such term leads to large degeneracy between cosmology and possible evolution 
I45I ) (Fig. El . 

UndoubtfuUy, despite its possible limitation, the determination of the Hubble diagram 
from SNIa heis led to a major and rapid change of paradigm in modern cosmology. 
However, this change has been possible because the previous situation was problematic. 
Although some observational indications were favoring a low density universe, the 
first dete ction s of fluctuations on degree scales were in conflict with open low density 
universe (1148). 



Fluctuations of the Cosmological background radiation The remarkable results of the 
WMAP experiment, with accurate measurements of the C; and additional measure- 
ments on the polarization, are often quoted as providing a direct evidence for an ac- 
celerating universe. This is incorrect: cosmological constraints established from CMB 
entirely rely on the spectrum shape assumptions, which is commonly assumed to be de- 
scribed by a single power law. Therefore the con clusi ons on the high precision obtained 
on cosmological parameters could be erroneous (|l25h . Indeed, relaxing this hypothesis, 
i.e. assuming a non power law power spectrum, it is possible to produce C; curve within 



23 




Fig. 5 Fitting the SNIa Hubble diagram with two free parameters, one being the cosmological 
constant in a flat cosmological model and the second being a parameter describing a possible 
time evolution of the luminosity of distant supernovae {Am{z) = K{tQ — t{z))) leads to the 
following constraints ll84f) . Contours are 1, 2 and 3 sigma regions. There is a severe degeneracy 
between the two parameters which prevents an unambiguous determination of the cosmological 
constant. A vanishing cosmological constant is entirely consistent with SNIa diagram provided 
that a significant but not unrealistic amount of evolution is allowed for. The major issue is that 
it is extremely difficult from a purely observational point of view to reject such a possibility. 



an Einstein de Sitter cosmological model which provides a fit as good as the concor- 
dance model. This is illustrated in figure [S] on which 3 models are compared to the 
WMAP data, two being Einstein de Sitter models. Such models not only reproduce 
the TT (temperature-temperature) spectrum, but are also extremely close in terms 
of ET (polarization-temperature) and EE (polarization-polarization) spectra. Further- 
more the matter power spectra are similar on scales probed by current galaxy surveys 
before the availability of the SDSS LRG sample. An un-clustered component of matter 
like a neutrino contribution or a quintessence field with ui ~ is necessary to obtain 
an acceptable amplitude of matter fluctuations on clusters scales ((33) . Such models re- 
quire a low Hubble constant ~ 46 km/s/Mpc. Such a value might be viewed as terribly 
at odd with canonical HST key program value ( ~ 72 km/s/Mpc) but is actually only 
^ 3(7 away from this value, this can certainly not be considered as a fatal problem for 
an Einstein-de Sitter universe. The introduction of a non-power law power spectrum 
might appear as unnatural. However, such a feat ure can be produced by some models 
of inflation in order to match the C; curve (jllj). Therefore the amplitude and shape 
of the CMB fluctuations as measured by WMAP is certainly a success for the A CDM 
model but cannot be regarded as a direct indication of the presence of dark energy. 



Large scale structure Once an Einstein de Sitter model is built in order to reproduce 
the CMB Ci , the amplitude of the matter fluctuations on large scales is set up and the 
measurement of the matter fluctuations on large scales in the present day universe is 
a critical way to distinguish models which are otherwise degenerated in their C;. The 
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Fig. 6 The TT spectrum of the first year WMAP data compared to three difi'erent models: 
one is the concordance, the two others are Einstein de Sitter models, one of which comprises 
neutrino contribution of ~ 10% corresponding to three degenerate families with m^ ~ 0.7eV. 
From Blanchard et al. 1133). 



comparison of the power spectrum from the SDSS LRG with the predicted spectra for 
Einstein de Sitter models is clearly in favor of the concordance model (|3J), see Fig. [T] 
One should add some caution here: it might be possible that the biasing mechanism 
leads to a power spectrum at small k (large scales) which is not proportional to the 
actual matter power spectrum (63), in which case the above comparison might not be 
a fatal failure of the Einstein de Sitter models. However, biasing mechanisms system- 
atically lead to a correlation function on large scales which is still proportional to the 
matter correlation function on large scales. Comparison of the correlation function on 
large scales is therefore less amb iguo us and its measurement should be unambiguously 
discriminant. Hunt and Sarkar l|ll2r ) have provided a comprehensive MCMC investi- 
gation of the Einstein de Sitter parameter space, finding models which acceptably fit 
the correlation function on scales below 70 h~ Mpc, but were nevertheless systemat- 
ically negative on scales of the BAO peak. This is a strong evidence that there is no 
way in an Einstein de Sitter universe to fit simultaneously the C; and the observed 
distribution of galaxies on large scales. This should be regarded as a remarkable suc- 
cess of the concordance cosmological model: although there were little doubts that this 
model could fit accurately most of the major existing observational facts in cosmology, 
the ability to produce predictions that are verified a posteriori is the signature of a 
satisfying scientific theory. Number counts of x-ray clusters were found to match the 
expectations of an EdS Universe and conflict with the concordance model. However, 
if the standard scaling of the M — T relation [3] is broken by some non-gravitational 
processes, the number counts and temperature distribution func tion could be altered 
in a way that the concordance model can accomodate data (|225r ). 



25 






10- 






10'" ^ 



10^ 



: ' 


' 


: 


: 




L . 


: 


U _j-^ '— ' 


■- -^ 


■ 


^ Jr^^ 


IcT^v 


■ 




'/ ^'^'rot^v^^ 


. 


..^'^ 


'^^tjbj^ 




y^ ■'V' 


^^*y^ 




y*^"*^ ■■'•' 


^Iml 


~ 


■ ■•^ 


•li^.. 


; 


.-■'y 


tWK 


■ 


.-■'y 


^MuK 


■ 


■■•■/' 


^!wl 


■ 


.■■x 


"^s 




..■V 


^ 


|k 


■.-■'^ 




tI^ 


^ 




^k ' 


■ 




\ ■ 


■ 




A 



0.001 



0.010 0.100 

k (h Mpc"') 



1.000 



Fig. 7 Data from the SDSS have allowed to measure the amplitude of galaxy fluctuations on 
large scales. In this respect, Luminous Red Galaxies (LRG) provided measurem ent of the power 
spectrum on the largest scales. Green crosses correspond to Tegmark et al. ||216| ) and black 
crosses correspond to the latest m easur ements of the power spectrum of LRG from the SDSS 
Data Release 5 by Percival et al. dlS^ ). The red continuous curve is the predicted spectrum 
for a typical concordance model, while the dotted and dashed lines correspond to the power 
spectrum for Ein stein de Sitter models consistent with the WMAP fluctuation angular power 
spectrum Ci l|3l : lll2f) . 



5.2 Einstein Cosmological constant or Vacuum contribution 

The most direct explanation one can provide to the cosmic acceleration it that it is 
due to a true cosmological constant appearing in the geometrical part of Einstein's 
equation, i.e. the left hand side of eq. U) However, it is much more common to believe 
that this term arises from some contribution to the energy-momentum tensor on the 
right hand side. As we have seen, from a classical point of view the vacuum might have 
a non zero density and behaves identically to a cosmological constant. In addition, 
quantum mechanics provides an intriguing hint in this direction. The possible energy 
levels of an harmonic oscillator are known to be: 



£^n = (n+ -)hv 



(67) 



and so the state of lowest possible energy, the zero point, is not zero but hv 12. This non- 
zero value is often noticed in standard text book of quantum mechanic, but, because 
observable quantities correspond to transitions from one state to an other one, is not 
regarded as being problematic. However, as soon as gravitational interaction has to be 
added, one cannot avoid to take the absolute energy into account. Summing all the 
contributions of modes of the fluctuations of the electromagnetic field up to some wave 
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number kc gives a density py. 



''" ATTk'^dk 1 , kc'^ 



The total contribution therefore depends on the cutoff scale kc- If the Planck scale is 
taken, this leads to a density which is something like 10 too large. For all energy 
scales in physics does py end up with an unacceptable large value, which looks like a 
fundamental problem [_|. An elegant solution to this problem is obtained from super- 
symmetry: the contribution to vacuum from fermion is negative and therefore with an 
equal number of modes in fermions and bosons one gets a cancellation. However, be- 
cause the supersymetry should be broken at energy below 1 TeV or so, the problem of 
the vacuum density is still not solved, even if its strength has been noticeabily reduced. 
An other proposition is that the vacuum actually behaves differently from l68t|40r ). 

Detecting the cosmological constant m a laboratory experiment? The above consider- 
ations might lead to a fascinating possibility. If the observed dark energy corresponds 
to an actual contribution of the electromagnetic field there will be a frequency fc as- 
sociated to the cut off in 1681 and the value associated with the present value of dark 
energy is: 

Vc ~ 1.7 lO^^Hz (69) 

Vacuum fluctu ation s produce noise which can be detected experimentally with Joseph- 
son junctions (|l28r ). it is therefore possible that the frequency cut-off kc could be 
meas ured in a laboratory experiment (jlSh . This point of view has been criticized by 
(|ll9r ). A further problem arises from the fact that the existence of a cut-off like 1691 
could mod ify the classical Casimir force(43), in a way which is excluded by existing 
data (|l58J). However, the origin of the noise current in the Josephson junction may 
have nothing to do with the electrom agne tic field in the device (JSOI) . Others arguments 
to reject this option have proposed ((ll85|)). 



5.3 The effect of inhomogeneities 

The standard formalism of Friedmann-Lemaitre models relies on the assumption that 
matter is uniformly distributed. Since the early work of Sachs and Wolfe (|200f l. per- 
turbations are commonly written as departure /i^^ from the Robertson- Walker metric 
element (7q,«: 

gal3 = 9al3+hal3- (70) 

This is a qualitative expression, as there are many ways to write perturbations depend- 
ing first on the choice of the initial form of the metric and coordinate systems. The 
analysis of linear perturbations (|/Iq,^| ^ (gaff)) is rich and not as trivial as one might 
naively think. However, the dynamics of relativistic linear perturbations to first order is 
well understood. Three distinct physical modes are possible: scalar perturbations which 
represent fluctuations of the density field p, vector perturbations which represent the 
vorticity of the velocity field and tensor perturbations which represent gravitational 
waves. Only this last term is obviously not present in Newtonian dynamics. This has 



^ It has been suggested that the zero-point fields should not be regarded as real, despite the 
fact that they are at the basis of the calculation of the Casimir effect {IGJ). 



27 

legitimated the use of Newtonian theory to describe evolution and average effect of 
perturbations on scales much smaller than the Hubble scale and for velocities much 
smaller than the speed of light. It is therefore common to use the Newtonian perturbed 
RW metric: 

ds'^ = -(1 + 2ilj)dt'^ + a{tf (1 - 2f/') -yijdx'dx^ (71) 

i/) being the Newtonian potential. However, the fact that we live in an inhomogeneous 
universe has always left open the possibility that average observable quantities may 
be definitively different from what can be obtained in a rigorous homogeneous world, 
this question di scuss ed in the context of relativistic world models, can be traced back 
to thirties (jG^; 12181 ). One early worry of this kind was about the average effect of 
gravitational lensing. It has been suggested that the average magnification from in- 
homogeneities in the Universe might produce an apparent redshift-magnitude relation 
which behaves differently from that in a homogeneous universe and would therefore 
bias inferred values of the cosmological parameters (JGy). S.Weinberg (2_3J) did provide 
a general argument to show that there is an integral constraint that guarantees the 
average relation in inhomogeneous worlds to be identical to the homogeneous case but 
that the way perturbations are handled may lead to lose this general integral con- 
straint on average quantities even to the first order in perturbation (while it is valid 
at all orders). Although Weinberg's theorem is extremely clever, it relies on some hy- 
potheses, one is that the universe still behaves as a FLRW model with its density being 
equa l to t he average density, and the subject has still been addressed in recent years 
(Im|;[i23;|ll;|2i). 

In this approach, the consequences of the presence of inhomogeneities were anal- 
ysed assuming the background evolution of the Universe was still following the usual 
Friedmann-Lemaitre equations. However, beyond the problem of properly evaluating 
the observable consequences of the presence of inhomogeneities (as illustrated by the 
calculation of the fiuctuations of the CMB), one fundamental question is whether the 
Friedmann-Lemaitre equations describing the expansion could be significantly altered 
by the presence of inhomogeneities, an effect named "back reaction" (JQJ). At first look 
one would naively expect that the average effect of linear perturbations to first order 
is zero and that to the second order one would get something like (for the case of 
perturbations around an Einstein de Sitter universe): 



(^ 



with : 



^-^{l + Fih^,)) (72) 



F(ft„^)(x {h^) (73) 



For the observable part of the universe, as we are within some perturbation with a 
typical size corresponding to the horizon, one expects that 

Pobs ~ P (l ± h^p) (74) 

(this is written in a very loose way that would probably horrify relativity experts, but 
it is done for the purpose of a qualitative illustration on how the problem is settled). 
Now astrophysical data show that (h) is always of the order of 10~ , or less, except in 
the vicinity of neutron stars and black holes, departures at this level would therefore 
not be very surprising. Actually the non linear collapse of structures like clusters or 
galaxies modestly enhances initial metric perturbations by a factor A ' , but lensing 
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measurements of galaxies and clusters directly prove that on these scales the metric 
perturbations remain tiny and one expects the linear approximation to be sufficient 
even when the structures are becoming non-linear (-^ 2> 1). Of course it is clear that 

some observables would be distorted at a level of /i ~ 10~ . However, the possible 
non negligible contribution of inhomogeneities to the global dynamics of the present 
universe, the so called "back reaction" has been proDosed I with the idea that the 
actual dynamics of the universe is much more distorted by non- linear terms than naively 
anticipated. Ironically enough, early propositions of back reaction were s uggested as a 
way to reduce or even cancel the presence of a large cosmological constant (|22ll : ll68l : l4ll ). 
while nowadays, the important question is whether this back reaction could explain the 
apparent cosmological constant. This perspective has attracted a lot of attention (see 
(|47|) for a recent review on this issue): it may even offer a possible solution to explain 
the acceleration of the universe, within the context of general relativity and without 
the introduction of any further exotic ingredient. T. Buchert (44) derived the exact 
dynamical equations governing the average size ajj of a domain D in an inhomogeneous 
matter distribution. These equations can be written in the following form : 

dx>\ _ 8nGp Q 71 , 



ax> J 3 6 6 

^)=-^ + | (76) 

aj) I 6 6 

for irrotational dust, and after dropping the cosmological constant for clarity. In this 
expression, 7?. is the average curvature over the domain and Q is an average quantity 
over the domain D: 

2=1 {{0^} - {ef - 2{a')) (77) 



being 3{ax>/aT>) a^nd a is the rate of shear a = 7^'^ij'^ij {'^ij being the shear tensor) 



2 1 

Given the analogy with the classical EFL equations, the possible consequence of 
back reaction has been discussed on this ground l|45l : 1461: 1471). 

Following Buchert's work several authors regard the back reaction problem as an 
essentially non pertu rbative question which cannot be solved without an appropriate 
treatment ((l9ll :ll3ll'). 

The question was handled from an other point of view by looking at se cond order 
effect from superhorizon perturbations generated during inflation ijlSl : llSOl ). The con- 
clusion was that some acceleration could be generated in this way. However, this result 
has been disproved (88; 98; 110; 99; 190). The same conclusion holds for superhorizon 
perturbations in presence of a scalar field l|l36l ). The possibility remains open that 
back reaction f rom subh orizon scales induces a modification of the equations governing 
the expansion (|l90l : |l3l|) and observational relations. However, from the astrophysical 
point of view what matters precisely is whether back reaction could alter significantly 
cosmological observations and fools the traditional interpretation. Clearly if the back 
reaction has to modify cosmological quantities by a factor h ~ 10^ the study of 
back reaction would remain of pure academic interest. Actually, modifications at the 
level of h ~ 10~ are expected from linear perturbations at the scale of our horizon. 
Indeed, (l23J) examined back reaction in a perturbative approach and concluded that 



^^ Tile idea that the small scale evolution may have a connection with large scale dynamics 
in relativistic cosmology can be traced back to Eddington in 1930 i6SD 
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gravitational energy of structure contributes as a back reaction term, but unless dark 
matter is essentially made of black holes, no appreciable change beyond the 10~ level 
is to be anticipated. Such a term should be called a weak back reaction, because it 
can be anticipated from simple dimensional arguments. The important question is : 
are there strong back reaction terms, i.e. terms which arise from non linear effects 
and which could contribute at a level much larger than the naive order of magnitude 
~ 10~ or even be of the order of being dynamically as important as the density p. At 
this level we are facing traditional difficulties when handling with relativistic effects. 
For instance, Buchert equation are written in a comoving frame in which particles are 
at rest. Therefore the non-zero "mass" associated with the kinetic energy of particles 
should appear as a "back reaction" term in the EFL equations. Identically the energy 
associated to a gravitational field has to be taken into account. These are weak back 
reaction terms. Using Newtonian perturbed RW metric (jlOJ) showed that back reac- 
tion term can entirely disappear from the second EFL equation [16] when the density is 
properly evaluated, but with a different choice of coordinates, and concluded that we 
are very close to a "no-go" theorem: "no cosmic acceleration occurs as a result of the 
nonlinear back reaction via averaging" (on this question a major source of complica- 
tion arises from the gauge issue and different authors may find apparently discrepant 
results because of the use of different gauges). The present formulation of Buchert 's 
equation is likely to depend explicitly on the choice of the coordinate system and there- 
fore the relevance to observations is unclear. It is not even clear whether a clever gauge 
transformation would not allow to suppress entirely this back reaction term. Similar 
conclusions were reached by (123}, who emphasized again that gravitational energy 



is a source of the gravity field which can appear as a "back reaction" ter m, i.e . is a 
weak back reaction term, and could be qualified as such, in agreement with (|234h . but 
which remains small in practice. In addition, they showed in some exact soluti ons t he 
absence of observable "back reaction". In this respect, the Swiss cheese-model (|206J) is 
an interesting exact solution in which the presence of inhomogeneities (that could even 
be non linear i.e. ft ~ 1) does not lead to any unexpected behavior i.e. in which there 
is no strong back reaction effect. From a second order perturbation analysis (|l23h also 
argued that as long the Newtonian conditions are satisfied: 

/ < -^ and u < c (78) 

"there always exist gauges where the metric differs from the RW metric only at or- 
der (ft)". Of course, when the size of the perturbations becomes comparable to the 
Hubble radius, one expects significant alteration of the standard relations (129;). A 
more troubling paper has been put on the ARXIV, but remains unpublished which : 
(jlTOl) investigated the solution of Buchert' equations in a exact spherical solution (the 
LTB that is described in the next section), actually a very interesting approach. They 
concluded that at some point "the un i verse starts accelerated expansion", a similar 
conclusion was drawn by (JS^; 1 15a : 1 16a : llSTI l , while it is clear that from spherical so- 
lution no physical acceleration is possible thanks to Birkhoff's theorem and as it has 
been confirmed by (jsl; \J7) ■ This is a direct indication that acceleration could appear 
in some quantities but which is by no means related to the observed acceleration. The 
question has fallen in the arena of subtle general relativity issues, but it seems that 
there are serious claims that back reaction cannot produce an acceleration compara- 
ble to what is requested from the observations. There are also convincing indications 
that no correction to standard EFL equations is expected above the amplitude of first 
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order perturbation 10~ . There is therefore an important "cosmological no-go theo- 
rem"which remains by now only a conjecture that for the standard CDM fluctuations 
spectrum, there is no signiflcant back reaction and the classical description is sufficient. 
It is by now the duty of advocates of real back reaction effete to demonstrate that in- 
teresting terms may exist beyond standard known linear and non-linear contributions. 
It remains important to examine whether weak back reaction could raise changes of 
the order of 10~ or more, as such an effect would start to be of the order of the 
precision anticipated for the determination of the cosmological parameters from future 
experiments. The fact that equation (|75|) looks like standard one , i.e. first EFL equa- 
tion, does not at all mean that terms appearing in it should be identified with similar 
looking standard cosmological parameters. These equations are valid for a given special 
choice of coordinates and it is therefore far from trivial to identify which observable 
consequences would follow. Even the term ^ should not a priori be identified to the 
standard Hubble expansion rate nor should ii-p be used in comparison with to —go • 
Strictly speaking when observers determine the Hubble constant Hq they compare the 
observed linear relation between luminosity distance and redshift from some galajcy 
sample, which are both observable quantities, and go is related to the leading second 
order term in z. Identically, when the cosmological density parameter (for instance) is 
determined by fitting the CMB Q , they do not proceed with the measurement of the 
actual density of the present day universe. A further problem is the fact that quantities 
entering eq.[75]and eq.[76]are averaged over some volume with some characteristic scale. 
They are therefore expected to vary with this scale. Indeed, virulent criticisms have 
been raised up against the (strong) back reaction program, pointing out that the metric 
eq. 1701 with Friedmann-Lemaitre equations have been fully successful to describe and 
to predict master pieces of observational cosmology, therefore any alternative model 
should demo nstra te its ability to reproduce these data as well with the same level of 
concordance ([117]). 

Claims have been made recent ly fo r the detection of observable ef fects tha t could 
be attributed to the back reaction (143). Some models have be en built |l6ll : l237l ) which 
were compared, successfully, to the main cosmological data l|l39l : llSSi ). in agreement 
with the "low" Hubble constant found by (214;). One can hope that this is the sign 
that the actual importance of back reaction will soon be clarified, although Cosmology 
is a fiel d we re there is a tradition for opponents to the standard paradigm not to resign 
easily (flTll ). 



5.4 Large scale void(s) 

After the discovery of an apparent acceleration from the SNIa Hubble diagram it has 
been emphasized that the Hubble diagram in its own was not a direct proof of an actual 
acceleration: even if the assumption that SNIa are standard candles is correct, an appar- 
ent acceleration could be obtained from the Hubble diagram within an inhomogeneous 
cosmological model fSO) when interpreted with the presumption of an homogeneous 
world. The possible existenc e of l arge scale structure up to 200 Mpc suggests that this 
option was to be considered (|219l ). The existence of very large scale structure has been 
a subject of regular claims in Cosmology (JGdulf). Such a possibility was regarded as a 
potential major failure of the standard vie w, and would even support the idea of a frac- 
tal distribution of matter on large scales l|93l : llSa : 1 1031 ) . The availability of very large 
scale surveys like the SDSS and the 2DF has almost entirely closed this issue: although 
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the visual impression from galaxy surveys might be that structures occasionally exist 
on very large scales, they are not the signature of average (r.m.s.) fluctuations on large 
scales greater than what could be anticipated for let say a standard ylCDM spectrum, 
when measured with appropriate tools like the correlation function or the power spec- 
trum. There is therefore no serious piece of data that suggests that the average level 
of fluctuations on large scales (as long as they remain smaller than the Hubble scale) 
is much larger than anticipated from the sta ndard pict ure. The anomalous character 
of some rare high fluctuations is still possible (|227l : 1112 ) and could for instance be due 
to some non- Gauss ian features of the primordial fluctuations, but the tendency of the 
human brain (|l8ll ) to identify insignificant patterns is a serious worry when dealing 
with these questions. 

Given that observations tell us that the universe is isotropic around us, if some 
large scale inhomogeneity has to exist and be in agreement with present day observa- 
tions it has to be nearly spherical and we should occupy a very specific place in this 
universe: for instance it may happen that we are close to the center of a near ly sp herical 
perturbation. Such a possibility might have some theoretical motivation ([1491). If we 
actually live inside a gigantic void, close to its center (to preserve the observed isotropy 
of the sky), this requires the abandonment of the Einstein cosmological principle that 
postulates a homogeneous matter distribution on large scales, and of the Copernican 
principle, that we are not in a special location in the universe. 

The homogeneity of the universe is directly observable in principle on scales much 
smaller than the Hubble radius, it is by no way obvious that it can actually be tested 
from observations on the largest scales (standard textbooks of relativity often men- 
tioned that it cannot be actually tested, but this is assuming a restricted set of ob- 
servable quantities). Detailed investigations of inhomogeneous models were carried out 
earlier, but not essentially with the purpose to provide an alternative explanation of 
the apparent acceleration, see for instance IjlSl ') for a description of various class of 
inhomogeneous cosmological solutions and references to earlier works on this subject. 

The exact solution of t he spheric al in homogeneity in general relativity was first 
published by Lemaitre (|l42l ). Tolman (|21g ) and Bondi (|42l ). aware of Lemaitre's work 
published studies on the same solutions. These solutions have been extensively used 
to examine closely whether the apparen t accelera tion could a ctual l y be due to such an 
inhomogeneous spherical solution (219l; ^^ [lli; [H; [13; 113; [lii; O)- More complex 
inhomogeneous world solutions to Einstein equ ations have also been proposed to repro- 
duce the observed SNIa Hubble diagram flla). Interestingly enough, (|27J) found that 
some "onion" structure could leave no strong apparent effect. Qualitatively speaking, 
a void is a region with a typical size jC where the density is lower than the average, 
in order for the initial singularity to happen at the same "initial time" in the past, 
the expansion rate in the region should be higher. This requirement is however not 
mandatory and it is possible to build isotropic inhomogeneous solutions with singular- 
ity starting at a different time (|49l '). Usually the "true" universe is an Einstein de Sitter 
Om = 1 a-nd the parameters inside the void are such that it looks like a low density 
universe universe. So that inside the void: 

^M ~ 0.25 and Hq ~ 0.72 (79) 

while asymptotically (outside the void): 

Hm = 1 and Ho ~ 0.50 (80) 
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the value of the Hubble constant being set by the "initial time" constraint. Roughly 
speaking in such models the distant universe has a Hubble constant lower than the local 
one, which is interpreted as an acceleration. The minimal size of the void to produce an 
apparent acceleration from supernovae at 2 ~ 0.5 is typically of the order 300h~ Mpc, 
or slightly less, i.e. up to 2 ~ 0.1 (j^). However, this kind of consideration is not sufficient 
to offer a satisfying alternative: the Hubble flow is known to be quite smooth and in 
agreement with an uniform expansion. A significant change of the Hubble constant over 
a volume of size smaller than the horizon size will lead to an appreciable change of 
the apparent Hubble constant. Actually, such effect has been found l|l20l ). who notice 
a Hubble constant 6.5% smaller on scales greater than 75h~ Mpc, a small number 
that might be due to a problem in the calibration of close supernovae (|56l). No other 
significant departure from an uniform Hubble constant has been recently reported, 
actually the Hubble diagram appears to be remark ably regu l ar fr om a few Mpc up to 



redshift where cosmological correction are needed (|l54l : Il33l : l214r i. Therefore, one can 
expect that voids witch size jC is significantly smaller than the Hubble radius will not 
mimic adequately the present day Hubble diagram. Appreciable difference would be 
easy to check in the near future. However, if the voids are large, comparable to the size 
of the Hubbl e rad ius, it is intuitive that a good match to the supernova data could 
be achieved. (|224r ) showed that in order to match closely the apparent acceleration of 
the Hubble diagram, a singularity at the origin should be present in the metric and 
other pathologies may exist. It is therefore far from being obvious that dust-filled LTB 
can reproduce the apparent acceleration in detail. Indeed, (55) concluded, that if one 
keeps the constraint that the curvature has to remain smooth around the origin, the 
Hubble diagram from an LTB will be different from the one in ylCDM and could be 
distinguished thereby with enough statistics at low redshit or intermediate redshift. 

Convincing examples of successive models have been obtained by (|7g;|95J), who also 
provided a public code : http://www.phys.au.dk/$\sim$haugboel/software.shtml 
For their best LTB model differences in the Hubble diagram between ylCDM and LTB 
appear to be less than 0.05m. Such differences are probably due to the choice of the 
analytical profile, and might therefore being reduced (or fitted to future data), even 
if this would be at the price of introducing further free parameters. The local Hubble 
constant is rather low {Hq ~ 60 km/s/Mpc) but still acceptable. It would be interest- 
ing to examine whether this behavior is generic or specific to their model. However, 
the confidence in the /ICDM relies on other pieces of evidence and the question of 
whether the CMB and the BAO feature can be reproduced has been also incorporated 
(|6|; y; [95I) . Because the outer cosmology is an Einstein de Sitter Universe with a small 
Hubble constant, an appropriate fit can be obtained relatively easily provided that the 
angular distance to the CMB is matched appropriately. These are rough evaluations, 
as a detailed formalism to deal with pertur bati ons in a LTB model is lacking. How- 
ever, this might soon become available (|54l : 12431 ). and it is reasonable to believe that 
an adequate fit can be obtained. It is nevertheless fair to say that these models are 
not yet as impressively good to fit cosmological data as the concordance is, perhaps 
because of their intrinsic complexity (but on the other hand, they contain many more 
degrees of freedom). The situation looks like that this type of hypotheses (large void) 
could be adjusted to the data and could not therefore be rejected fro m ob servations if 
enough tuning is allowed. However, fortunately, this is not the case! |223|) have shown 
that in a LTB model, the time drift of the cosmological redshift can be different from 
what it is in homogeneous worlds, providing a possible test of the Copernican principle. 
The kinematic Sunyaev-Zel'dovich (kSZ) effect on distant clusters will be significantly 
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different for distant clusters offering a different method of testing LTB (|96l ). Anotlier 
impressive tool for constraining LTB models has been proposed by (|48r ): we know that 
the intergalactic medium which probably contains most of the baryon is highly ionized 
up to redshifts greater than 5; CMB photons are scattered by the electrons of this 
plasma, so photons we collect from the CMB are a mixture of photons which have 
travelled straight towards us and of photons which have been scattered at lower red- 
shifts, if these electrons are within the void region, they are scattering a very different 
CMB because they are moving rapidly. In addition, even electrons which are out of the 
void will see a different CMB sky because the void itself and will be a further source 
of a distorted CMB. Therefore contrary to the homogeneous case, an observer in the 
center will observe a combination of black-bodies with different temperatures, resulting 
in a distorted spectrum. Present day limits on possible spectral distortion of the CMB 
already provide interesting constraints and could become critically more stringent if 
limits could be improved by an order of magnitude (|48l l . 



5.5 Quintessence 

A true cosmological constant is an object that is largely unwanted from the theoretical 
point of view. Ratra and Peebles introduced the idea that the acceleration could be 
due to the presence of some scalar field dominating present day density of the Universe 
(1192; 1182 ). Quintessence names a scalar field that coupled to gravity. The canonical 
Lagrangian of a scalar field is : 

C=l{d,^f~VW (81) 

the first term being the kinetic energy {X = ^{di(f>) ), and the second one the potential. 
The stress-energy tensor has a form identical to that of an ideal fluid with pressure 
and density given by : 

p = C = ^{d,^f^V(^) (82) 

p = 2XdxC-£ = ^{d,4>f + V{<t}) (83) 

For a field which is spatially homogeneous the equation of state parameter w is: 

i0^-T/(,^) 

w = - — -=—. 84) 

P i02+T/(^) 

which remains > — 1. The equation driving the field in a FL world is : 



, + 3H^^-^{^) (85) 



The original scenario was proposed l|l92f) with a potential of the form: 



4-(-Q 



M' 
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Let us suppose that we are at some late time when radiation or matter dominates. Then 

1 2 
the expansion factor is a{t) oc f" with n = —, —. A power law solution of equation H85II 

cj) cc t^ behaves as: 

cj> ex t^ (87) 

and the ratio of the density of the scalar field to the total density is: 



— i- oc a"(°+2) 
P 

So that if Q > the density of the scalar field will be dominant. Another important 
property is that solution [57] is an attractor in that it describes the late solution of a 
large class of initial conditions. If the kinetic term becomes small at a later epoch in 
1841 the field behaves like a cosmological constant. This does not necessary solve the 
coincidence problem, but at least alleviates it. 

Dark energy has become the subject of intense theoretical efforts, which lead to 
investigate different forms of the Lagrangian. One first generalisation consists in mod- 
ifying the potential V((j)). A more radical approach is to modify the kinetic term. In 
fc— essence models IjlOtllH). inspired from fc— inflation ([9|), the Lagrangian is written in 
a general form : 

C = p{^,X) (89) 

and contains therefore a non canonical kinetic term. 

Few constraints can be added. It can be requested that the density p remains pos- 
itive, although there is no decisive argument against a negative cosmological constant 
and so this constraint might not be necessary. On the other hand such a term would 
lead to a rapid big crunch, which is unwanted. Another possible constraint is that the 
square of the sound's speed remains positive (9): 

2 _ 9X£ _ Of (QQ. 

d d^ 

(,X stands for -—rp and , ATX „ ). 

Such models easily lead to large variations of Cs (while quintessence based on (|8ip 
automatically leads to Cs — 1) and one can have Cs > 1, this is not problematic because 
the theory remains Lorentz invariant, so there is no violation of causality (1781 : Il2rl I . 
Within such models the equation of state parameter w could take any value ((|l63)) 
and so one can have: 

w < -1. (91) 

Such dark energy is referred to as phantom (ghost) energy. A simple example is ob- 
tained by changing the sign of the kinetic energy term in l81l 

So from the astrophysical point of view w can be regarded as a free function which 
values should be determined from observations. From the fundamental point of view 
V{4)) or the Lagrangian itself are the quantities to be determined and therefore observa- 
tions have to be used to constrain directly these quantities, rather than the parameter 
w{z). Finally, let us mention that more details on scalar field models of dark energy 
can be found in the recent review (158 



For a different point of view see (|37 
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5.6 New gravity law on large scales 

The existence of a large scale acceleration in the Universe is a serious indication that 
our present knowledge of gravity is actually incomplete. Although the introduction of 
a new component, quintessence, might provide the explanation for the acceleration, 
i.e. a new term in the energy momentum tensor, the historical cosmological constant 
introduced a new term in the other side of Einstein equation. Therefore it might well 
be that the geometrical term is more complex that we used to believe. This option of 
modifying gravity to account for the acceleration is briefly introduced in this section. 
More details can be found in recent reviews on this topic (|62l : ll72r) . 

While General Relativity can be derived from the Einstein-Hilbert action: 

5 = / dx^R (92) 

where R is the Ricci scalar, one can try a more general action, leading to the so called 
f{R) theories: 

dx'^fiR) (93) 

where / is an unknown function. Such theories have to pass solar system tests as well 
as pulsar chronometry tests on small scales and cosmological tests on large scales. In 
practice such theories met considerable problems, and it is far from being clear that 
viable models can be built in this way. An alternative procedure, is the so called Palatini 
formalism, in which the variation principle is modified: the metric and the connection 
are regarded as independent quantities, and the resulting equations are in general 
different but for f{R) oc R. Within this formalism, (7) have examined constraints on a 
specific model: 

f[R) = Rll + a(^\ j (94) 

where the second term is essentially a correction term to the classical Einstein-Hilbert 
Lagrangian. Acceptable regions of the parameter space encompass the standard ylCDM 
model, but these models differ from standard cosmological models, as the relations 

between J?^,, Qm, qo and the growth rate are different. 

a In a 
A more radical way to modify gravity in our world is through the idea of higher 

dimensi onal space. Such propositions can be traced back to 1919 with the work of 

Kaluza (|l22r ) who proposed an unification scheme for gravity and electromagnetism 

within a fifth dimensional space. (|l27t ) pointed out the interest of having a fifth compact 

dimension to avoid observational constraints on large additional dimension. Superstring 

theory and supergravity theories possess remarkable properties in higher dimensional 

space and have therefore deserved strong attention from theorists. Modern versions are 

known as braneworlds or brane cosmology (8; 24). Our 3+1 world lies on the brane, 

while the remaining space is the bulk. Matter, with pressure and density, is present 

only in the brane, but gravity is present in all dimensions. The vacuum energy in the 

brane provides a tension term a and there is an other vacuum energy in the bulk 

Aq- The general action now contains terms involving an equivalent of the Ricci scalar 

corresponding to gravity in higher dimension. The inferred equations describing the 
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expansion on the brane, the generalized FL equations, contain more new terms which 
could be non- linear and which are related to the properties of the bulk \2a ): 

the third additional term behaves like radiation but comes from the bulk and is some- 
ti me ca lled the dark radiation on which BBN provides stringent constraints (|24l : l25l : 
SEH). 

Direct modification of the FL equations have also been proposed through Cardas- 
sian models (|90|): 

H^ ^Ap + Bp" (96) 

which does not contain a vacuum contribution, although in models with n = the 
second term of the right hand side behaves identically to a yl term. Such a form for the 
EFL equation can arise from some particular fluid properties of dark matter or from 
higher dimensions (|l02l ). 

In DGP models, extra dimensions are infinite (i641 ), and the effective action contains 
explicitly a 4-D Einstein-Hilbert term i.e. on the brane in addition to the 5-D term on 
the bulk. This introduces a scale re and two distinct regimes appear: on scales smaller 
than re gravity essentially results from the 4-D term and classical GR is recovered 
(although high precision tests might lead to some differences (115g: il'i11 '): on larger scales 
the gravity is "leaking" the expansion is eventually accelerated. The FL equation is 
replaced by ((|6l|)): 

H^^(jP+11 + e^y (97) 

VVSMJ,, 4r2 2rcJ 

with e — ±1. From this it is clear that at early times =- ^ — =■ allows to recover 

3M|,; 4r2 

the classical EFL regime, while at late times the accelerated expansion is recovered 

(provided that e = +1). Detailed predictions of DGP models may not be identical 

to thos e of y lCD M, and recent investigations found some tension between data and 

theory (|l99l : 11981 : 1821 ) . although the predictions are not as straightforward as in the 

standard model. From the astrophysical point of view, an interesting aspect of these 

classes of models in which gravity is modified is that relations between cosmological 

parameters are not identical to t hose in the ylCDM and the evolution of the growth 

factor is expected to be different (| 1461 1 . with a possible dependence on scales, allowing 

for possible discrimination between various possible origins for the acceleration (12221 : 

H). 

A more radical option has emerged in recent years : given the fact that the ac- 
celerated expansion clearly needs some revision of our knowledge of the gravitational 
sector, could it be that the problem of dark matter itself is due to the break of stan- 
dard gravity laws at finite distance? This possibility has been introduced in an empiric 
way by Milgr om in order to explain the rotation curves of galaxies without invoking 
dark matter (165). A fully relativistic theory has been built by Bekenstein (19j which 
behaves like Milgrom' law in the weak field regime: the TeVeS (tensor-vector-theory) 
theory. Quite remarkably, it has been shown recently that this theory can reproduce 
the observed cosmic acceleration, large scale power spectrum (at the time of the work) 
and CMB with comparable success as ylCDM L209, ). The present status of this theoy is 
not as satisfactory as with the standard concordance model : this theory did not lead to 
specific predictions which have been verified a posteriori, and it needs the introduction 
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of a new type of fields in piiysics in the form of vector fields. However, MOND is a clear 
example of an alternative view, which differs drastically in its fundamental ingredients 
and illustrates the fact that our understanding of the gravitation sector relevant to 
cosmology might be more limited than commonly assumed. 

6 The area of precision Cosmology 



Parameter 


Vanilla 


Vanilla + n^ 


Vanilla + w 


Vanilla + Ok + w 


nth-^ 


0.0227 ± 0.0005 


0.0227 ± 0.0006 


0.0228 ± 0.0006 


0.0227 ±0.0005 


Hch^ 


0.112 ±0.003 


0.109 ±0.005 


0.109 ±0.005 


0.109 ± 0.005 


e 


1.042 ± 0.003 


1.042 ±0.003 


1.042 ± 0.003 


1.042 ± 0.003 


T 


0.085 ±0.017 


0.088 ±0.017 


0.087 ±0.017 


0.088 ± 0.017 


Us 


0.963 ±0.012 


0.964 ±0.013 


0.967 ±0.014 


0.964 ±0.014 


n. 





-0.005 ± 0.007 





-0.005 ± 0.0121 


w 


-1 


-1 


-0.965 ± 0.056 


-1.003 ±0.102 


^x 


0.738 ±0.015 


0.735 ±0.016 


0.739 ±0.014 


0.733 ± 0.020 


Age 


13.7 ±0.1 


13.9 ±0.4 


13.7 ±0.1 


13.9 ±0.6 


Om 


0.262 ±0.015 


0.270 ±0.019 


0.261 ±0.020 


0.272 ± 0.029 


0-8 


0.806 ± 0.023 


0.791 ±0.030 


0.816 ±0.014 


0.788 ± 0.042 


Zre 


10.9 ±1.4 


11.0±1.5 


11.0 ±1.5 


11.0 ±1.4 


h 


0.716 ±0.014 


0.699 ±0.028 


0.713 ±0.015 


0.698 ±0.037 



Table 1 Summary of the mean values and 68% confidence intervals for the parameters con- 
strained from CMB, SNIa and BAO for different models {9 is the ratio of sound horizon to 
angular diameter distance). These constraints are quite tight, most of them are below 5%, and 
are stable when additional degrees of freedom are added to the model {w.fifS), adapted from 



After the discovery of the fluctuations of the cosmological background by COBE the 
perspective to achieve precision measurements of the angular power spectrum of these 
fluctuations appears to be within reach and two satellite experiments were designed to 
reach this target: WMAP and Planck. This has opened a new avenue for Cosmology 
to beneflt from high precision measurements with well controlled systematics. Indeed, 
WMAP delivered data of high accuracy, allowing for high precision estimations of the 
cosmological parameters (|212l : 1631. Furthermore, the consistency of SNIa, BAO an d 
CMB data allows reliable accurate estimations by combining constraints (|l32l : Il33l i. 
The possibility of reaching high precision measurements with different techniques (dis- 
tant SNIa with SNAP, properties of distant x-ray Clusters as proposed by Panoramrx 
some years ago or the recent WFXT ([1691), or with the SZ sample of clusters expected 
with Planck and other CMB experiment like SPT, full sky weak lensing surveys with 
DUNE, or the more recently proposed redshift surveys of hundred millions of galaxies 
with SPACE or ADEPT) has not only reinforced the perspective to determine cos- 
mological parameters with high precision, but will also allow to investigate the very 
nature of dark energy. The Dark Energy Task Force report (2!) presented a summary 
of the different levels of progress expected for various projects. Ultimate constraints 
will be obtained by combination of various experiments. It is interesting to summa- 
rize what present day data allow. There are many possible data to combine, but the 
three most currently used are SNIa, CMB and measure of the large scale distribution, 
either through the power spectrum or through the correlation function. A combina- 
tion of these three data sets leads to tight constraints on the minimal A Cold Dark 
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Matter model. These constraints are summarized in Table |6] adapted from (|8J). Dif- 
ferent groups rec ently got simi lar results with approaches which slightly differ in the 
technical details l| 133 : 1 1331 : 12051 ), so these numbers can be regarded as very robust. As 
one can see, most of the cosmological parameters describing our Friedmann-Lemaitre 
universe are constrained to an accuracy better than 5%. Furthermore, when the pa- 
rameter space is enlarged the constraints remain essentially unchanged. This calls for 
some caution. Quoted uncertainties reflects statistical uncertainties. Unidentified sys- 
tematics are the critical issues in this topic. Variations in published values of ag from 
various approaches (CMB, Clusters, weak lensing) have provided an illustration that 
systematics uncertainties could alter preferred values beyond statistical uncertainties. 
On the other hand, given possible systematics which have been identified until now, 
it is likely that future estimations will remain within the two sigma domains. First, 
when combining only two probes, one already gets tight constraints which are within 
this range. A second argument is that when one allows for additional parameters (free 
equation of state parameters w, curvature, non-zero neutrino mass, tensor contribution, 
...) preferred values and interval ranges are not changed by much. These are indica- 
tions that we are already in the precision area of cosmology: present day estimations 
of cosmologcail parameters are likely not to change by much in the near future and 
investigations of the nature of dark energy will need extremely accurate control of sys- 
tematics. Whether the necessary investments will be valuable for the astronomical field 
is a subject of debate l|236l l. 



7 Conclusions 

The Copernician model of the world was the first revolution of a series in the con- 
struction of modern cosmology, and the discovery of the accelerated expansion being 
the latest in date. Theoretical considerations have always been a source of remarkable 
observational investigation and Cosmology has always benefited from the confronta- 
tion of models with observations. Since the thirties, the big bang picture, the modern 
version of Lemaitre's primeval atom has been remarkably successful, based on simple 
assumptions and physics laws that have been validated by accurate experimental re- 
sults. Although alternative theories have been developed, these alternative were based 
on hypothetical unknown physics advocated to interpret cosmological observations. 
None of these alternative theories has produced predictions that have been comforted 
a posteriori. Rather new observations in agreement with predictions of the big bang 
picture necessitated deep revision of the unorthodox view, at the cost of rather ad 
hoc assumptions added to fit the new observations. The situation has evolved when 
the standard picture has necessitated the introduction of new ingredients, first dark 
matter and more recently dark energy. The very nature of these new ingredients, which 
are supposed to dominate the mean density of the universe has not been established 
by direct laboratory experiments, nor by astronomical observations, and this situation 
may some time lead to the question whether cosmologists have not introduced new 
aethers. We had the opportunity to see that the situation is not so. The introduc- 
tion of -cold- non-baryonic dark matter has led to specific predictions, the amplitude 
of the fluctuations of the cosmological background on various angular scales, which 
were verified with high accuracy precision. The presence of dark energy has lead to a 
specific prediction, the shape of the power spectrum on large scales, which hcis been 
verified a posteriori. Although, the inclusion of a cosmological constant was concomi- 
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tant to general relativity, the actual origin of dark energy remains totally unknown and 
the presence of dark energy in the present day universe represents probably the most 
fundamental and unexpected new element in modern physics. 
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